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We consider the Lemaitre-Tolman-Bondi metric with an inhomogeneous viscous fluid source sat- 
isfying the equation of state of an interactive mixture of radiation and matter. Assuming conditions 
prior to the decouphng era, we apply Extended Irreversible Thermodynamcs (EIT) to this mix- 
ture. Using the full transport equation of EIT we show that the relaxation time of shear viscosity 
can be several orders of magnitude larger than the Thomson collision time between photons and 
electrons. A comparison with the "truncated" transport equation for these models reveals that the 
latter cannot describe properly the decoupling of matter and radiation. 
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H ' The relaxation time r for the decay of any dissipative flux in a hydrodynamic system is often misunderstood and 

1 I , uncritically assumed to vanish (an unphysical assumption) or be of the order of the relevant collision time tcoi of the 
' particles making up the fluid -for a recent review see . The root of this confusion may be traced to the great success 

of the Fourier-Navier-Stokes theory of dissipative transport This venerable theory is very useful for engineering 
applications, but it ignores altogether the relaxation times taken by the system to return to thermodynamic equilibrium 

I state once the gradients of temperature, fluid velocity and so on are suddenly turned off, and therefore its range of 

I . validity is limited (at most) to steady-state situations Therefore, this theory displays a serious drawback: the 

^ ' prediction of instantaneous propagation of the dissipative signals which, on physical grounds, should be limited by the 

, speed of sound. This fact becomes important when dealing with relativistic self-gravitating fluid distributions such 

J^D ' as a collapsing star, as well as for the early Universe (see [Q, ||] and references therein). Obviously r and tcoi ought 

, to be related to one another, however there is no generic rule on how to find such a relationship since it depends in 

■ each case on the particular system under consideration. 

y—( I In a previous paper we studied the matter-radiation decoupling for an expanding inhomogeneous universe described 

. by new exact solutions ||^ characterized by the Lemaitre-Tolman-Bondi (LTB) metric with flat space sections [^|J|] 

' and a viscous fluid source. The dissipative quantities involved (the shear stresses) in these models are well behaved 

qh, and comply with relativistic causality and stability in the sense of Israel and Stewart [p[-pd| . As stressed in Q , the 

■ study of inhomogeneous cosmologies is a well motivated and justified endeavor, especially in view of the results of 
Mustapha et al. who show that inhomogeneous but spherically symmetric models can be compatible with current GBR 

• observations This paper aims at showing that the solution found in ||] implies that the relaxation time of the 

. 5^ I dissipative shear stress of the cosmic matter-radiation mixture prior decoupling is several orders of magnitud larger 

' than the collision time between photons and electrons. In section II we present our model, recall the exact solution 

?H derived in , and compare the aforesaid relaxation time with the collision time between photons and electrons before 

5^ the matter-radiation decoupling is achieved. Section III summarizes our conclusions. 
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II. THE LTB DISSIPATIVE MODEL 



Cosmic matter sources in the period from the end of cosmic nucleosynthesis to the decoupling of matter and radiation 
(roughly the temperature range lO^K < T < lO^K) can be characterized as a mixture of ultra-relativistic and non- 
relativistic particles ("radiation" and "matter"), subjected to a tight interaction via various radiative processes (hence 
the name "radiative" era). The usual approach for studying this era is to consider a FLRW space-time background 
where matter sources are described by: (a) equilibrium kinetic theory, (b) gauge invariant perturbations, or (c) by 
hydrodynamical models, which, in general, fail to incorporate a physically plausible description of the matter-radiation 
interaction. Examination of these matter sources in an inhomogeneous hydrodynamic context requires an imperfect 
fluid whose state variables are compatible with the equation of state of a mixture of a non-relativistic gas (matter) 
and a extreme relativistic one (radiation). In Q we derived a class of exact solutions along these lines, so that: (a) 
matter-radiation radiative interaction was modeled by a dissipative shear-stress tensor, (b) only the evolution along 
the radiative era was considered, (c) the background space-time was described by the LTB metric with flat spacelike 
slices: 
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ds^ = -c^dt^ + dr^ + [de^ + sm^{0)d(f>^] , Y = Y{t, r). (1) 
For a comoving 4- velocity u°- — cSf in (|l|), the expansion scalar and shear tensor are 

Y' 2Y 1 f Y Y'\ 

<a^0 = ^ + — , fT%-diag[0,-2a,a,a], a ^ - iy - - \ (2) 

where a prime denotes derivative with respect to r and Y = u°'Y^a = Y^t- 
The stress-energy tensor of the matter-radiation mixture is 

T"^ = pu'^v!' + ph"^ + n'^'' [h"^ = c-^u^u^ + „ = 0) , (3) 

where the most general form for the shear-stress for (1) is given by: diag [0, — 2P, P, P], with P = P{t,r) to be 
determined by the field equations. The relation between p and p follows by assuming that the matter component 
is "dust" (the low temperature limit of a classical ideal gas, hence its pressure can be neglected) and therefore the 
hydrostatic pressure is furnished by the radiation only, i.e. 

where the superscripts (m) and (r) denote the particle number densities of "matter" and "radiation" . The shear 
dissipative pressure tensor is governed by the transport equation of EIT 



1 / r , 
—Tri — — u 

2 '\Tij 



4 , . 

+ 2?7 aab = 0, with 77 = r;,^^, = -p^'^V , (4) 



where r is the relaxation time for the shear-stress and rj is the coefficient of shear viscosity. The subscript rg in 
T] indicates that we use the form of this coefficient provided by kinetic theory for the electron-photon interaction 
(the "radiative gas" model). This transport equation is compatible with relativistic causality and stability, and is 
supported by kinetic theory, statistical fluctuation theory, and information theory ||3|,P - 111 . In its turn the expression 
for rj is well-known in the literature and it can be derived from different standpoints |3j^3[. The entropy per particle 
takes the form 

ni 

where a = —T/{2ri^^^^ ) — — 5/(8p^''-'), s'-'^-' is the equilibrium photon entropy per baryon, and we have used the fact that 
the number of particles of each fluid entering the mixture is independently conserved: (n'^™^M°);a — {n'-^^u°')-a ~ 0, 
as well as n^™^ jn'^'^^ <^ 1. The positive-definite nature of s and r above is crucial for the theory to comply with 
relativistic causality, notwithstanding the unjustified and widespread belief that r vanishes altogether or is at most 
of the order of the collision time. 

A. Exact Solution and Density Contrasts 

Integration of Einstein's field equations (see details in [^) leads to the following exact solution 

^VM(^-t^) = Vy+^(?/-2e)-VTT^(l-2e), (5) 
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where the definitions 



were used along with 



_ kM _W _Y 

^ = W' = M' Y 



M = / fJr'Y^Y^dr, p<r^ EE mc'n^r' , f^^Y^Y^dr, p« ^ 3n^' k^T, , 
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so that p'l^\ pf^ define the initial densities of the non-relativistic and relativistic components of the mixture, re- 
spectively. The subscript "i" refers to some initial hypersurface to be specified below, and k, denotes Einstein's 
gravitational constant. 

For the inhomogeneous solutions derived above it is convenient to characterize initial conditions by introducing 



suitably averaged initial densities: 



(r) 
Pi 



and density contrast parameters Aj , Aj defined by 



Am) 



1 + A, 



1 + A 



as well as the ancillary functions 



Y'/Y 



I* = 1 



(1-r) 

3(1 + Ar)^ 



$ = 1 



(1 - 4r) 

3(1 + Ar')' 



that appear in the fmictional forms of the state variables: n'™) = n'f"'^ / (y^T) , n'''^ — ii'p /{y'^T), T = Ti'^/y and 
P = {l/2)n[^^kgTi^/{y^T). These functions are linked to the density contrasts by 



r = 1 + 3AA|"' + SBAf ' , 



(6) 



* = 1 



AA 



BA 



(r) 



(7) 



-AA^r + (1 - a: 



l + A 

where the quantities A and B are known functions of y 



(t) 



(8) 



A = 



1 
3"? 



= (1 - 4e - 8e2) - (y^ _ _ ge^) 



B = - 
V 



(l + 2e)^^-(y + 26) 
V 1 + e 



It should be noted that the solutions must satisfy the restriction F > so that the particle number densities do not 
become negative and that no shell-crossing singularities occur. 

At this point it is expedient to point out first that no expression exists that gives t solely in terms of the macroscopic 
variables (p, p, T, n, etc.), i.e., there is not such thing as an "equation of state" for r. Secondly, using the kinetic 
theory of gases (in general) the collision time (not r) could in principle be written as a collision integral, but that 
expression would be mathematically untractable. This is why we resort to (Q) to get r in terms of the quantities 
introduced above 



a 



1 



171 p2 
256 



(9) 



This expression is justified as long as it behaves as a relaxation parameter for the interactive matter-radiation mixture 
in the theoretical framework of KIT. 

Given a set of initial conditions specified by e, A^ , A^ , on some initial hypersurface, the forms of the state and 
geometric variables as functions of y and the chosen initial conditions are fully determined. However, for the solutions 
to be physically meaningful they must comply with the following set of physical restrictions. 



16t 



>0, 



(10) 



(consistent with the condition that the entropy production be positive). 
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* > 0, 



cr$ < 0, 



(11) 



to ensure that r, p and T are all positive. Likewise the concavity and stability of s demand that 



f > 0, 



S _ 2crr \P^ 



> 



1 + 




(12) 



The collision time {t^ ) for the Thomson scattering of photons by electrons (the dominant process of shear transport 
in the radiative era) is given from the Saha equation by 




where ct^ , Bq, rrie and h are, respectively, the Thomson scattering cross section, the hydrogen atom binding energy, the 
electron mass and Planck's constant. For details of the derivation of (13) see Before matter-radiation decoupling 

< tn (where tn = 3/0 denotes the expansion time, i.e., the "Hubble time"). For t-y > tn matter and radiation no 
longer interact between each other. 

Choosing as initial time t = ti roughly corresponding to T = 10^ Kelvin, the evolution of the models crucially 
depends on initial conditions specified by the useful quantity A-**^ = | A-'^^ — A^™^ (formally analogous to the definition 
of adiabatic perturbations). As shown in (see also Figure 3 of this paper), the collision time tj is quite insensitive to 

(s) ..... . 

the values of A^ given by initial conditions, therefore for some some y = ?/d we find that overcomes tn, i.e. there 
is always a decoupling hypersurface. Taking the standard decoupling temperature Tq w 4 x 10"^, numerical evaluation 
of Td = (IO^/j/d) 5'(j/_D, A-'*-', A-'^-') where 4' is given by (|^) leads to yo ~ 10^-^. Also, observational constraints 
on the decoupHng hypersurface [STeq/Teq]D ~ 10~^ further restrict the maximum values of both A,-™^ and A-'^^ to 

about 10~^. However, r is very sensitive to A^*'', for example, only if A^"^ = (i.e. for adiabatic initial conditions) 
T has an adequate asymptotic evolution, in the sense that as y — s- cx) (and so, t — s- oo) it grows monotonically (as 
J- —j. y3/2-j g^^^ g ^ in this limit, though r in this case remains smaller than Ih for all the evolution as ?/ — s- cx), a very 
different behavior from that of tj which overtakes Ih at the decoupling surface. The relaxation time for adiabatic 

( s) 

initial conditions is then unphysical near the decoupling surface. For non-adiabatic intial conditions, A^ ^ 0, if 
A|[j'A|^'' > and (T<f> < hold, we can have cr —> for non-vanishing therefore r diverges and overtakes tn, a 

behavior that is qualitatively analogous to that of t^. In order to have t overtaking tn at ?/o w lO^ "*, we find that A^^^ 
cannot be greater than approximately 10~^ (quasi-adiabatic initial conditions). Therefore, under the assumptions of 
EIT and using the full transport equation (Q) and r given by (^), we find that initial conditions do exist so that t 
has a physicaly reasonable form near the decoupling surface. Figures 2, 3 and 5 illustrate this conclussion. 



On the initial hypersurfce = 5*^ = 1, $i « A,;*^^ and Ti is evaluated to be Ti « 10^ seconds. Likewise [t^]i « 2227 
seconds. As a consequence [r/tj^Ji ~ 5/16 while [tj/tH]i ~ 10"^. As the decoupling hypersurface is approached 
(i.e., y yn) ^ and F slightly depart from unity as depicted in Figure 1. On that hypersurface one has tj{yD) « 

Tivo) ~ 10^'^ seconds. Figure 2 shows how sensitive is the latter quantity to We see therefore that contrary 

to the usual belief, the model examined in ||^ shows that t is much larger than for most of the radiative era and 
the timescales r, tj become comparable only as the decoupling hypersurface is approached. This behavior is shown 

in Figure 3 depicting the ratios T/tn and t^/tu as well their dependence on |A|*''|. 



In many cases, just for the sake of mathematical simplicity, the transport equation (Q) governing the evolution of 
the dissipative stress tensor has been replaced in by the more manageable Maxwell-Cattaneo-like expression, namely 
the so called "truncated" transport equation 



B. Comparing the different times 



C. The truncated transport equation 



Tiled Khi + Uab + 277 Cab 



0. 



(14) 
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an equation still complying with relativistic causality and stability in the sense mentioned earlier, and justified when 
the four-dimensional gradient of the combination tu'^ /{Trj) becomes negligibly small. The models presented here 
provide a theoretical context in which to test how suitable can this truncated transport equation be for studyin g th e 
evolution of an interactive matter-radiation mixture. To do this consider the relaxation time that follows from 



P 4 + SA^''^ - 4r 



P + Ipa 4m + 3A|"^ - 4r] - i(T[31(4 + SA,^"^) - 15r] ' 



(15) 



where Tmc is the relaxation time associated with the Maxwell-Cattaneo equation. Comparing (^) with the relaxation 
time (|) that emerges from the full theory (i.e. from (M)) 

tA I qa('') Ar\fA I qa('') -r^ 

(16) 



P _ 1 (4 + 3Af^ -4r)(4 + 3Af' -r) 



P+lpa- ^P 4a 4[4^3^M_4pj2^ mp2 

we can see that r = Tmc(l + Tmc{p/p))^^ , therefore since p/p < (as long as ^ > 0, which is always the case) we 
have r > Tmc- Now two possibilities arise: 

1. Adiabatic condition A|'*'' = 0, in this case Figure 4 shows that t,„c < tn for the whole evolution. 

2. Quasi-adiabatic condition A^*'' ^ 0, in this case as seen from Figure 5, Tmc does not diverge nor overtakes tn 
at any time in the whole evolution, while r diverges and overtakes tn for y > yo ^ 10^'^. 

This together with the fact that close to the decoupling surface becomes larger than Tmc it follows that the 
Maxwell-Cattaneo equation (|l^) applied to the LTB models derived here and under the constraints imposed on the 

initial density contrasts (A^-'^'' « 10^^ and 10~^ < |A|'*''| < 10^^) by the observed CMB anisotropics, does not provide 
a relaxation time that behaves reasonably in conditions near the matter-radiation decoupling. 

It is widely assumed nowadays that the bulk of the matter of the Universe (about ninety percent) is not baryonic 
but belongs to some hitherto unknown form of "dark matter" (see e.g. |l^) -though voices of dissent can be heard [ p^ . 
We have considered a baryon-photon mixture and so have not taken this into account in our calculations. However, 
we can roughly verify to what extent our results would change had we incorporated some form of this missing matter. 
If we assume that it is mostly non-relativistic and that it amounts to 20 times the baryon abundance, the inclusion of 
such missing matter corresponds to decreasing the value of the parameter e by a factor of 20 (from 10"^ down to 50). 
As shown in Figure 6, carrying on this modification does not imply a dramatical change in our results. Further, it is 
worthwhile remarking that the introduction of a small cosmological constant or a quintessence fluid would not make 
sense for the models we are dealing with. Since these forms of "dark energy" contributions (as it is nowadays in vogue) 
are used to explain the alleged accelerated expansion of the present day Universe, they should be sub-dominants at 
the epcoch contemplated in this paper (see e.g. |p^). 

III. CONCLUDING REMARKS 

It is worthwhile mentioning that the main results of this paper: (a) inadequacy of the truncated equation ( p^ to 
provide a relaxation time that behaves reasonably near matter-radiation decoupling, and (b) that r is much larger 
than t.y for most of the radiative era, strictly apply for the specific model examined in based on the LTB metric. 
These results also rely on the expression of r obtained in terms of the transport equations (|^) and (p^. However, 
we argue that in spite of these limitations these results can complement and enrich the existing discussion in the 
literature concerning the adequacy of truncated transport equations of dissipative stresses in a cosmological context 
-see [0 and references therein. Our analysis suggests that, in any case, these equations should be used with due 
care. We believe it is important to stress that (as far as we are aware) these important features of the equations 
of irreversible thermodynamics have not been examined previously in the context of inhomogeneous metrics able to 
accommodate dissipative shear stresses. 
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LIST OF CAPTION FOR FIGURES 



Figure 1 

The plot depicts the functions ^ and F, given by (6) and (7), in terms of log;Lo(2/) ^^'^ logiod^i*^ I)' from the initial 
hypersurface y = 1 corresponding to Tj = 10^ Kelvin. Notice how these functions are almost unity for all the evolution 
range. 

Figure 2 

Logarithmic plot of r vs [A^*-*! for y — y^-^ corresponding to T = T^. Note that r diverges for |A^*^| w 10~^. Since 
f-y also diverges as T — > T^, the behavior of r for this value of |A^*'| is reasonable. 

Figure 3 

Logarithmic plot depicting the ratios Tjt^ (shaded surface) and t-^jij^ (non-shaded surface) vs the parameters logif){y) 
and lA^"-*!. The former ratio is far larger than the latter for most of the evolution range. For |£>eZtof ^1 w IQ-^ the 
two ratios become comparable in magnitude near the decoupling surface {t-y/t^ = 1). 
Figure 4 

Logarithmic plots of the ratios t/Ih (top surface) and Tmc/tjj (bottom surface) vs the parameters y and |A^''-'| = 

(4/3) I A,^"^ I under adiabatic conditions (lA,^**^] = 0). Neither one of these relaxation times diverges nor overtakes Ih, 
therefore they are not physically realistic near the decoupling surface tj/tn = 1- 

Figure 5 

The same plot as figure 4 but under quasi-adiabatic conditions. The ratio r/tn diverges for logioy « 2.4 while the 
the ratio r^c/tH remains finite for all y, therefore the full transport equation yields a reasonable relaxation time while 
the truncated equation does not. 

Figure 6 

This figure compares r/t^ (A), Tmc/t^ (B) and t-y/t^ (C) vs \ogiQ{y). The thick curves were obtained assuming 

T = 10^ Kelvin and an initial photon to baryon ratio jrif^^ = 10^, hence e « P^l^ I Pi^^ = 10'^. If we assume that 
CDM does not interact with the photon gas and that it constitutes most of non-relativistic matter, say 20 times the 
abundance of baryons, then we would need to use e = 50, lower by a factor of 20. The thin curves were obtained using 
this value of e with |A-*^| = 10~^ (upper curves) and \l^f'\ = 10~^ (lower curves). Notice that these two curves are 
indistinguishable in (C), thus indicating that tj is rather insensitive to changes in jA^^'*!. The changes in e due to the 
inclusion of the rest-mass of CDM are also minimal, shifting the best value of |A^*'''| to an order of magnitude less. 
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